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Abstract 

We give a rational expression for the subresultants of n + 1 generic polynomials 
fl, ■ ■ ■ , fn+i m n variables as a function of the coordinates of the common roots 
of fi, . . . , f n and their evaluation in / n +i- We present a simple technique to prove 
our results, giving new proofs and generalizing the classical Poisson product for- 
mula for the projective resultant, as well as the expressions of Hong for univariate 
subresultants in roots. 
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1 Introduction 

The classical Poisson product formula for resultants of univariate polynomials 
can be stated as follows: if / and g are two univariate polynomials of degrees 
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d\ and d<i respectively, with g = bd 2 (x — £1) • • • (x — £d 2 ), then the resultant of 
/ and g can be expressed as 

Res(f,g) = (-l) d ^b%Y[f(Q. (1) 

3=1 

The main result of this paper is a generalization of Formula (1) for univariate 
and multivariate subresultants (see Theorems 2.2 and 3.2). Although most 
of the results in the univariate case already appeared in [17,18,23,9], here we 
present simple techniques that enable us to reobtain them (see Theorem 2.2 
and Corollary 2.6) and allow us to generalize them to the multivariate case. 

Resultants and subresultants of two univariate polynomials go back to Leibniz, 
Euler, Bezout and Jacobi. We refer to [12] for historical references. In their 
modern form, subresultants were introduced by Sylvester in [26]. They have 
been used to give an efficient and parallelizable algorithm for computing the 
greatest common divisor of two polynomials [7,3,15,10,19,20,25]. More recently 
they were also applied in symbolic-numeric computation [11,31,22,30]. 

Multivariate resultants were mainly introduced by Macaulay in [24], after 
earlier work by Euler, Sylvester and Cayley, while multivariate subresultants 
were first defined by Gonzalez- Vega in [13,14], generalizing Habicht's method 
[16]. The notion of subresultants that we use in the present paper was intro- 
duced by Chardin in [5]. It works as follows: let . . . , be a system of 
generic homogeneous polynomials in K[ ] of degrees di = deg(/f) 

with parametric coefficients, where s < n + 1 and K is the coefficient field 
of fiT--,fs - Let 7i<zi,...,d s : N — > N be the Hilbert function of a complete 
intersection given by s homogeneous polynomials in n + 1 variables of de- 
grees di,...,d s . Fix t G N and let S be a set of Ti.d lt ...,d s (t) monomials of 
degree t. The subresultant As is a polynomial in K whose degree in the co- 
efficients of is TC dl ,...,d i _ u d i+1 ,...,d 3 (t - d^ for i = 1, . . . , s, having the follow- 
ing universal property: A s vanishes at a particular coefficient specialization 
fi, fg £ C[x , . . . , x n ] if and only if I t U <S does not generate the space of 
all forms of degree t. Here, I t is the degree t part of the ideal generated by the 
#'s (see [5]). 

The constructions in [13,5] generalize the classical univariate subresultants 
in the sense that they provide the coefficients of certain polynomials in I t , 
which in the univariate case include the greatest common divisor of two given 
polynomials. 

Theoretical properties and applications of multivariate subresultants are active 
areas of research. A series of recent publications explored: their application to 
solve zero dimensional [14] and over-constrained polynomial systems [28], in 
the inverse parametrization problem of rational surfaces [1]; their irreducibility 
and connection with residual resultants [2] ; the generalization of their universal 
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properties to the affine well-constrained case [8]; as well as generalizations of 
matrix constructions for subresultants [27]. 

Multivariate subresultants also encapsulate as a particular case the classical 
projective resultant Res(fi,...,f% +1 ), which is defined to be an irreducible 
polynomial in the coefficients of the f^s which vanishes at a particular coeffi- 
cient specialization /f, f% +1 e C[xq, . . . , x n ] if and only if f% +1 have 
a common root in the complex projective space Pg. 

There is an affine interpretation of the resultant that can be stated as follows: 
Set 

fi //*(!, a^i, • • -, x n), fi ■= f t h (0,xi,.. .,x n ), i = l,...n + 1. 

Due to Bezout's Theorem, the cardinality of the set 

V(fi, ...,/„): {iiT': f,(0 = / 2 (0 = . . . = f n (0 = 0} 

equals di...d n (here, overline denotes algebraic closure), and the classical 
Poisson product formula [29,6,21], which generalizes (1), states that the fol- 
lowing identity holds in K 

Res(/f , f h n+1 ) = Res(/i, . . . JJ** 1 II WO- (2) 

Cev(/i,...,/„) 

In order to make this formula a generalization of (1), we have to define 
resultants of non-homogeneous polynomials. The obvious generalization is 
Res(<7i, . . . , g n +i) '■= Res(<7ij ■■■i9%+i)i where is the homogenization of gj. 
The same extension to affine polynomials holds for subresultants. It should 
also be mentioned that the Poisson formula (2) is a particular case of the 
determinant of a multiplication map in a quotient ring (see [21, Prop. 2.7]). 

In Theorem 3.2 we generalize (2) and give an expression for any multivariate 
subresultant as a ratio of two determinants times a function of the coefficients 
of f\, . . . , f n . The determinant in the denominator is a Vandermonde type de- 
terminant depending on the common roots of /i, . . . , f n , while the determinant 
in the numerator depends on evaluations of the common roots of /i, in 
the last polynomial f n +i- 

The paper is structured as follows: in Section 2, we present in detail the uni- 
variate case, showing how to derive with our techniques Hong's expressions for 
subresultants of two univariate polynomials in the roots of one of them and 
the coefficients of the other. The details in the univariate case are essential for 
the generalization to the multivariate case: they allow to identify the extra- 
neous factor which is non-trivial in the multivariate case and they also allow 
to handle the generality of the monomial sets appearing in the definition of 
multivariate subresultants. In Section 3, we deal with the general case. 
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In order to keep coherence with the classical literature and previous works, the 
presentation in the univariate case is done in the traditional way, i.e. for non- 
homogeneous polynomials, while in the multivariate case the reader should be 
aware that the notions involve homogeneous polynomials. 

Acknowledgements: This work started during a visit of the first two authors 
to North Carolina State University in the Fall Semester 2004. We would like 
to thank the Department of Mathematics for their hospitality and for the 
stimulating working atmosphere it provided. A previous version of this paper 
was presented in Mega'05. We thank the anonymous referees of that previous 
version for introducing us to discrete Wronskians as presented in [23] and for 
pointing us other related references. 



2 The univariate case 



Classical scalar and polynomial subresultants 

We review here the definition and some well-known properties of the classical 
univariate resultant and scalar subresultants and polynomial subresultants. 

Let / = a dl x dl + . . . + a and g := b d . 2 x d2 + . . . + b = b d2 (x - £i) • • • (x - £ d2 ), 
a dl 7^ 0, bd 2 7^ 0, be two polynomials of degrees d\ and di respectively with 
coefficients in a field K and roots in the algebraic closure K. 

The scalar subresultant of / and g is defined for < j < j < min{<ii, g^} 
as the following determinant: 



:= dct 



di + d 2 — 2k 

0>di ' ' ' " ' Ofe+l-(d2-fe-l) a j-(d 2 -k-l) 



do — k 



(3) 



>d 2 



Wl-(di-fc-l) fy-(di-fc-l) 



di — k 



M 2 



Jk+i 



where ag = b e = for £ < 0. 
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The subresultant polynomial SreSfc(/, (?) is defined for < k < min{di, d 2 } as 

Sres k (f,g) :=£Sj?V. 



When A; = 0, Sreso(/, g) = coincides with the classical resultant Kes(f,g) 
which arose historically when checking if / and g have a common factor: 



gcd(/, g) = 1 



Res(/,«?)^0. 



In an analogous way, the scalar subresultants satisfy the following property: 



deg gcd(/, g) = k 



Sf ] = for < £ < k and S k K> + 0, 



(k) 



and the polynomial subresultants sresk(f,g) are determinant expressions for 
modified remainders in the Euclidean algorithm. In particular, for the first k 
such that ^ 0, the monic gcd of / and g satisfies: 

gcd(/,^) = (5f ) )- 1 sres fc (/^). 



There is a generalization of the univariate Poisson product formula (1) for the 
polynomial subresultant sres fc (/, g), as shown by Hong in [17, Th. 3.1], see also 
[23, Formula 9.3.2] and [9, Sec. 5]: 



STes k (f,g) = (-l)^- k ^-% d d l- k 







&m) • • 


■ (x-uxv 


tdi— k— 1 
SI 




£d 2 —k- 
' Srf 2 


7(u) 




e? • 


to 

• Sd 2 






e?- 1 • 


td 2 -\ 
■ Sd 2 





(4) 



(Here the sign is due to the fact that we consider / on the roots of g instead 
of g on the roots of / as done in [17].) 
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Notations: 



As we mentioned earlier, most of the results we obtain in this section are not 
new. However, we consider important to illustrate our technique by applying 
it to the univariate case, since it helps to understand its generalization to the 
multivariate setting. The choices of notations we made here are accordingly 
motivated by their coherence with the multivariate case. They correspond to 
Chardin's notion of subresultants [5] applied to the univariate slight 
generalization of the usual notion of scalar subresultants. 

• / := a + a\X + . . . + a dl x dl and g := b + b±x + . . . + bd 2 x d2 in K[x], where 
K := Q(a , . . . , a dl , b 0: . . . , bd 2 ), with a , . . . , a di: b 0: . . . , bd 2 algebraically in- 
dependent variables over Q (representing the indeterminate coefficients of 
two generic polynomials / and g of degrees d\ and d 2 respectively). 

• • • • ) £,d 2 } denotes the set of roots of g in K (recall that overline denotes 
algebraic closure), and Vd 2 '■= det(£j _1 )i<jj<<2 2 the Vandermonde determi- 
nant associated to this set. 

• For any j G Z, K[x]j : = {0} U {/ G K[x\ : deg/ < j}. Note that if j < 0, 
then K[x]j = {0}. 

• We set t G Z such that 0<t<d 1 + d 2 — 1, and let t* := max{ci 2 — l,t}. 

• M f G ijr(*-<*i+i)x(f+i) and M g G K(*-^+ 1 )x( t *+ 1 ) denote the transposes of 
the matrices in the monomial bases of the composition of the Sylvester 
multiplication maps and the inclusion K[x] t — > K[x] t *: 



[i f : K[x\ t - dl -> K[x\r fi g : K[x\ t _ d2 -> 

and 

r p ^ x l3 g(x) 



x 



x a f(x) 



x h 



where the monomials indexing the rows and columns of these matrices are 
ordered "increasingly" 1, x, x 2 , . . . . Namely 



Mf 



ao - ■ .a<n 



ao • • -a dl 







bo- ■ -b d2 


















b ...b d . 2 





Note that if t < d\ then Mf = (the empty matrix), and if t < d 2 then 
M g = 0. 
We set 



k := t + 1 - dim (K[x] t - dl ^ - dim {K[x\ t - d ^j 

= t + 1 - max{0, t - di + 1} - max{0, t - d 2 + 1} 
= t + 1 - max{0, i - di + 1} - (r - d 2 + 1). 



(5) 



Note that k > since t < c?i + c? 2 — 1. 

5 := {a; 71 , . . . , x lk ; < 71 < • • ■ < 7^ < t} C K[x] t , a fixed set of k 
monomials of degree bounded by t. 
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• sg(«S) := (— l) a where a is a number of transpositions needed to bring 
(1, X, x 2 , ... , x 1 *) to 

^71 ™7fc ~ t+1 t** 1 t ^.71-1 ^.71+1 ^.72-1 „72+l 

I «X • • • ^ JU j JU 2***7 *-*--' 5 5 5***5 5 5 *** 5 5 5 *** 5 / * 

• As : = Ag\f, g) denotes the order t subresultant of f,g with respect to S, 
i.e. the determinant of the matrix whose max{0,£ — d\ + 1} first rows are 
Mf, whose max{0,t — d 2 + 1} following rows are M g and from which one 
deletes the k + t* — t columns indexed by S U {x t+1 , ...,#**}. 

Remark 2.1 The order t subresultant of f,g with respect to S coincides (up 
to a sign) with the scalar subresultant when making special choices of t and 
S: 

(1) When * = d x + d 2 - 1, then fc = t + 1 - d 2 - di = and S = 0. In 
that case, from the definitions of Res(/, g) and A® one gets that A© = 
(-l)**Res(/,<7). 

(2) For < k < min{di,d 2 } and t := d\ + d 2 — k — 1, we can take Sj := 
{x l ,0 < i < k,i 7^ j}. In that case, from the definition of A$. and (3) 
one gets that A s . = (_i)(di-k)(«fe-*) 1 S'CO. 



The main statement of this section corresponds to (a slight generalization 
of) Hong's theorem [18, Th. 3.1]. It expresses A$ as the ratio of discrete 
Wronskians: we refer to [23, Sec. 9.3] for an introduction to the subject. Here 
we present a new simple proof of this result, that we generalize in the next 
section to the multivariate setting. 

Theorem 2.2 Let f,gE K[x] and . . . , £d 2 } be the set of roots of g in K . 
Then, under the previous notations, for any fixed t, < t < d\ + d 2 — 1, and 



for any S = {x 



71 



x 



} C K[x] t of cardinality k, with k defined in (5), the 



order t subresultant As of f, g with respect to S satisfies: 



sg(S) 6; 



where 



i-yi 



t*-d 2 +l 

d 2 



£71 
Sd 2 



\Os\ 




£<2 2 1 f{£>d 2 ) 



e K 



■d 2 X d 2 
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Proof. First, Os is a square matrix since by (5) we have 



d 2 = k + (t* - t) + max{0, t-di + 1}. 



Let Is £ j ff( fe + t * _t ) x (**+!) be the transpose of the matrix of the immersion 
of the ^-vector space generated by S U {x t+1 , . . . ,x f *} into if[x] t * (Is is an 
identity (k + t* — t)-square matrix plugged into (t* + 1) zero columns), and 
set 



M s :-- 



(6) 



Since it is straightforward to check by (5) that we have 

k + t* — t + max{0, t-d 1 + l} + max{0, t - d 2 + 1} = t* + 1, 
therefore Ms is a (t* + l)-square matrix. 



Furthermore, it is immediate to verify that \Ms\ 
to prove that \M S \ = b d \~ d2+l \O s \/V d2 . 



sg(S)A 5 , and we are left 



We set 

to to 
si ' ' ' Sd 2 

and we observe that Vd 2 = \Vd 2 \- Now, we perform the product M s V^ 2 : 



w (f+l)xd 2 T/ . 














-^(t*+i)x(r+i) 




Id 









* 









r 


* 










* 




Id 







6d 2 o 




* bd 2 



□ 



Therefore |M 5 | V d2 = b%- d2+1 \O s \, which proves the Theorem. 
The following examples illustrate how the formula works in a couple of cases. 



Example 2.3 d\ — 5, d 2 — 2, t — 4. Now we have t — t*, k = 2, and 



M / = 0, M g 



b b 1 b 2 
b 6i 6 2 
6 6i &2 



Set 5 := {rr, a; 4 }. Here As does not coincide with any of the scalar sub- 
resultants S 2 \ < j < 2. However, it is straightforward to check that 
As = b bf — 6q6 2 - On the other hand, since sg(«S) = 1, by Theorem 2.2 we 
have that 



sg(5) b\ 





6 


6 


2+1 


i\ 


6 4 




l 


l 




6 


6 



^ 66 



£f £ 3 
; 6-6 



= fr!66[(6 + 6) 2 -66] 

= ( 6o / &2 ) [( 6l / 6a )2 _ (6o/62 )]. 



□ 



Next example deals with a case when t < d 2 in which case we need to use 
t* — d 2 — 1 instead of t. 



Example 2.4 di = 2, d 2 = 5, t = 3. Here 
associated to this value of k are S^ 2 " 1 ' ; c ' ( 
for t = 3 < d 2 we have t* = d 2 — 1 



2. The scalar subresultants 

a 2) S^ 1 ' = a 2°i an( i *S , 2°' ) = a 2°0) while 
4. Thus we have 



do fli 02 

a ai a 2 



M o = 0. 



For 5 := A s = a|, and Theorem 2.2 still works in this case: since 

sg(«S) = 1 and &5 _5+1 = 1, one has 



1 1 


l 


1 1 




1 


l l 


l 


l 




2 6 6 6 




6 


6 6 


6 


6 


£i e 


4 £4 £-4 £4 
2 S3 S4 S5 




6 4 


£4 t4 
S2 S3 


6 4 


£4 
S5 


7(6) 7(6) /(6) /(6) /(6) 




6 2 


£2 £2 
S2 S3 


a 


6 2 


6/(6) 6/(6) 6/(6) 6/(6) 6/(6) 


- = a 2 r 


6 3 


6 3 6 3 


a 


6 3 




11111 






l 


l l 


i 


l 




6 6 6 6 6 




6 


6 6 


6 


6 




£•2 £2 £2 £2 £2 
Si S2 S3 S4 S5 




6 2 


S2 S3 


6 2 


6 2 




£3 £3 £3 £3 £3 
SI S2 S3 S4 S5 




6 3 


6 6 


6 3 


6 3 




£4 £4 £4 £4 £4 

Si S2 S3 S4 S5 






6 6 


6 4 


6 4 
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□ 



We end this section by showing how simple it is to derive from Theorem 2.2 
both the Poisson product formula (1) and Hong's formula (4) for subresultant 
polynomials in roots, together with its generalization to a larger class of deter- 
minant polynomials that we call here generalized subresultant polynomials. 

Observation 2.5 (Poisson product formula) Applying the previous theorem 
to Remark 2.1(1), one obtains 

Res(/,s) = (-l) dld2 A 



f-1 



\did 2 idi 
' °d 2 





' ^d 2 f(^d 2 ) 






• drV(a) 




& ■ 








g-1 . 


td 2 -\ 





d 2 



f -iy id2 b% n KM 
3=1 



□ 



Observation 2.6 ([17, Th. 3.1], [23, Id. 9.3.2]) We derive Hong's Formula 
(4) applying Theorem 2.2 to Remark 2.1(2): 



sres fc (/, g) = ^ S t 



3=0 



= (_l)(di-fc)(d2-fc) ^ a s .x j 

3=0 

(-^^^^^^r^Esg^oi^ 

3=0 



X J . 



We observe that in this case t* = t, sg(Sk) = sg{l, . . . , x k x } = 1 and sg(«S,-) 
( — l) fc ~ J , and thus, by column expansion of the determinant we get: 



£sg(<S;)|(9 5j 
3=0 



X/ 



(-1)* 

{-l) k x 

(—l) k x k 




e } 



ei* _ *"7(ei) 



to 

Sd 2 
£<2 2 

tfe 
W 2 

£d 2 f(£d 2/ 



£dl k 1 f(€d 2 ) 
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= (-1)* 



1 

tf-ztf 



tf-xtt 1 



rf 2 
fc-1 



£rf 2 ^2 



dr fe_1 /(u) 



td^—k—X 
?rf 2 



□ 



One can straightforwardly generalize Hong's result to a larger class of deter- 
minant polynomials 

k 



s(x) :=£A 5j 



X 



1i 



(7) 



3=0 



corresponding to an arbitrary set of monomials S := {x lj , < j < k} C if[x]t 
and Sj := S \ {x 7j } , where c?2 < £ < ^i + ^2 — 1 and k : = c?2 — max{0, t — c?i + 1 } . 
We call such a polynomial a generalized subresultant polynomial. 

The usual proof that shows that sresfc(/, g) belongs to the ideal (/, g) generated 
by / and g extends to showing that s G (/, g) and the following expression in 
terms of roots holds (we omit the proof which is essentially the same than the 
proof of Observation 2.6). 

Corollary 2.7 Let f,g e K[x] and s(x) be the generalized subresultant poly- 
nomial defined in (7). Then, we have 



(x 7!-70_ ^-70^70 
c 7fc"7fc-i _£7fc-7fc-l\£7fc-l 

^ 2 " fe_1 /(6) 



( x 7fe-7fe-i _^fc _ Tfc-l)^fc-l 
?rf 2 /(^d 2 ) 

£^2 f(€ch) 



□ 
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3 The multivariate case 

In this section we generalize Theorem 2.2 to Chardin's multivariate subresul- 
tants [5], after introducing the notations we need. 

Notations: 

• For n G N and 1 < % < n + 1, 

\a\<di 

where a = (a±, . . . , a n ) G (Z> ) n , x a := a;" 1 • • • x" n , |a| — a± + • • • + a n , and 
-K" := Q(aj Q , 1 < i < n + 1, | or j < di), with a ia algebraically independent 
variables over Q (representing the indeterminate coefficients of n + 1 generic 
polynomials in n variables fi of degrees di respectively). 

• For any j G Z, K[x]j := K[x u 4 = {0} U {/ G K[x] : deg / < j}. 

• We set t G N, p := (d 1 — 1) H h (d„ - 1) and £* := max{p, t}. 

• k :— T~td 1 ...d n+1 (t), the Hilbert function at £ of a regular sequence of n + 1 
homogeneous polynomials in n + 1 variables of degrees d±, . . . , d n+i , i.e. 

k := #{x a : |a| < t, at < di, 1 < i < n, andt — \a\ < d n+ i}. 

• S := {x 71 , . . . , x 7fe } C K\x\ t a se t of /c monomials of degree bounded by t. 

• For 1 < % < n + 1 , 

7£j := {x a , |a| < t — di, a 3 - < dj for j < i}. 
We observe that for 1 < % < n, 

#(7?-j) = #{x a , | ok | < t, otj < dj for j < i and aij > di}, 

and 

#(^ n+ i) = #{x Q , |a| < i, «j < dj- Vj and i - \a\ > d n+1 }. 
Therefore 

iV:= + = dim K K[x] t = k + #(K 1 ) + --- + #(K n+1 ). (8) 

• In particular, we denote 1Z n+ i =: {x^ 1 , . . . ,y^ r }, where r := #(1Z n+ i) and 
we observe that 

k + r = #{x a , \a\ < t, a, < dj Vj} = dim K[x\ t /{h, ...,/„) D K[x] t . (9) 

• For j > 0, Tj := Tid 1 ...d n (j), the Hilbert function at j of a regular sequence 
of n homogeneous polynomials in n variables of degrees d\, . . . , d n , i.e. 

Tj := #{x a : \a\ = j, < di for 1 < % < n}. 
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We note that tj — if j > p. 
For j > 0, 




{x Q : \a\ = j, di < di for 1 < i < n} for < j < t — d n+ \ + 1. 



any set of Tj monomials of degree j for j > max{0, t — d n+ i + 1}, 



(10) 



See Remark 3.3 for a discussion on the definition of Tj. 
T := Uj> Tj and T* := Uf =t+1 Tj. We note that #T = d, where d := 
di • • • d n is the Bezout number, the number of common solutions of fi, . . . , f n 
in K n , and that T* = if t* = t, i.e. if t > p. 

In particular, we denote T = {x ai , . . . ,x ad }, and we assume that T* = 
{x% . . . ,x as }, the first s := #(T*) elements of T. 

if [x] tj * denotes the if -vector space generated by K [x] t U T* and A^* : = 
dim(if[x] ti *). 

For 1 < i < n + 1, Mf. £ x dlTa( - n ^ xN * denotes the transpose of the matrix in 
the monomial bases of the composition between the Sylvester multiplication 
map and the inclusion K [x] 4 — > K [x] tj #: 



For later convenience we order the monomial basis of if [x] tj * in such a way 

that all monomials in T precede the monomials in if [x] tj * \ T. 

M5 e f('( 7V - fc ) x ( Af - fc ) denotes the Macaulay-Chardin matrix obtained from 



by deleting the columns indexed by the monomials in S U T*. 
Following [24,5], we define the extraneous factor £ (t) as the determinant of 
the square sub matrix of (11) whose rows are indexed by all those monomials 
x a G TZi, 1 < % < n, such that t — dj — | ct | > d n+ i or there exists j > % with 
ctj > dj, and whose columns are indexed by those x a such that t — \a\ > d n+ i 
and for some index i, a.i > d^ or such that there exist at least two different 
indexes 1 < i, j < n with ctj > di, otj > dj. It is straightforward to verify 
that this is really a square matrix. An important property of £(t) is that it 
neither depends on the coefficients of f n+ \ nor on S. 

A s := Ag\(fi, . . . , fn + \) denotes the order t subresultant of . . . , f£ +1 



with respect to S h := {x 71 ^^ 1 ', . . . ,x 7fc a;^^ fc '}. Here, fj 1 denotes the ho- 
mogenization of fi by the variable x n+1 . It turns out that by [4] we have 



x a 1— > X Q /j 



(11) 



M 



"n+1 
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± 



I Ms I 
£{t) 



;i2) 



• For 1 < i < n, f i is the homogeneous component of degree di of fi, and 
A-Tj := A^f ) (f 1 , . . . , f n ) is the order j subresultant of f\, . . . , f n with respect 

to ' ' _ f) 

• {Ci) • • • ! £d} denotes the set of all common roots of /i, . . . , f n in K , and 
Vr := det(^" i )i<jj<d, the generalized Vandermonde determinant associated 
to T. 

Remark 3.1 The order t subresultant given in (12) generalizes both the uni- 
variate case and the usual multivariate projective resultant as defined for in- 
stance in [6, Th. 2.3]. 

(1) When n = 1 and t < d± + d% — 1, there are no rows and columns of (11) 
satisfying the condition that contributes to the extraneous factor £ t , and 
thus £(t) = 1. Therefore of (12) coincides with the univariate order t 
subresultant of / and g with respect to S defined in Section 2. 

(2) When t > p + d n+ ±, then k = since ci\ < d±, . . . , a n < d n imply \a\ < p, 
thus t — \a\ > d n+ \. Therefore S := 0. In that case we recover Macaulay's 
construction [24, Th. p.9 and Th. 4] and A 5 = ±Res(/ 1 A , . . . , f r h l+1 ). 

We are ready now to state the main result of the paper, the multivariate 
generalization of Theorem 2.2. 

Theorem 3.2 Let f\, . . . , f n+1 e K[x\ and . . . , £ d } be the set of common 
roots of fi, . . . , f n in TC 1 . Then, under the previous notations, for any t e Z>o 
and for any S = {x 71 , . . . ,x 7fc } C K[x\ t °f cardinality k = Tid 1 ...d n+1 (t), the 
order t subresultant A$ satisfies: 



= ± I II 

\j=t-d n+ i+l 




(13) 



where 



o. 



a 1 ■ 


£71 
Sd 




Clk 
Sd 


er • 


Sd 


SI 


Sd 


CfVn+l^l) - 


' CdVn+l^d) 


_4 r f r /n+l(6) ' 


• £d7n+l(£d) 



dxd 
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Proof. First we check that Os is a square matrix, i.e. that d = k + s + r. This 
is clear by Formula (9) since 

d-s = #(T) - #(T*) = #(T x T*) = #{x a , |a| < f , a 3 < d j Vj} = k + r. 
In this proof the monomial basis {x^ 1 , . . . , x 5Ar * } of if[x] t) * is ordered such 
as was specified in the notations (monomials in T precede the rest of the 
monomials in K[x] t ^). 

Like in the univariate case, we define Is £ K^ k+s ' )xN * as the transpose of the 
matrix of the immersion of the if-vector space generated by SUT* into if [x] ti * 
in the monomial bases. We set 



Ms :-- 




G K 



N* xN* 



(Ms is a square matrix by (8) and since N* = N + dim(T*).) 

Furthermore, it is immediate to verify that \M S \ = ± \M S \ = ±£(t) A 5 , where 
£(t) denotes the extraneous factor that has been introduced in (12). 

We set 



Vm* 



- a 

■ ■ Sd 



e K , and V d := 








Vn* 






Id 



e K 



N*xN* 



and we observe that Vr = |Vd|- We perform the product M s V d . 



M s V d 











id 



^3 


* 


^3 


* 




M h 









M' f 

Jn 


£j*fn+l{£j) 


* 
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where M 



i ._ 



M' h 



Jn 



is the submatrix of 



M fl 



M 



fn 



with the same number 



of rows and whose columns are indexed by all monomials in x a e if [x] t) * \T = 
K[x\ t \ (T \ T*) = K[x\ t \ T . It is immediate to verify that M' is a square 
matrix since, again by (9), #(7^i)+- • •+#(72 n ) = N-k-r = iV-#(T \T*) = 
iV* - d. 



We recall that #(T \T*) = #{x a , |a| < t, on < d { Vz}, and therefore M' is the 
Macaulay-Chardin matrix associated to the computation of A^ r , (fi, . . . , /„), 
the order t subresultant of fi, . . . , f n with respect to T \ T*. 

To conclude the proof we are left to prove that 



\m'\ = ±e(t) i n ^ 

\j=t-d n+ i+l 



This was proven in [24, p. 14] (see also the proof of [4, Lem. 1] and [8, Thm. 5.2]). 
For the reader's convenience, we rewrite the proof here. 

We reorganize the matrix M' as follows: we recall that the columns correspond 
to monomials x a G K[x\ t \T and we index the columns by graded descending 
order, first all monomials of degree t in if [x] t \T, then all monomials of degree 
t — 1 in K [x] t \ T, and so on, up to all monomials of degree t — d n+ i + 1. 
Finally, we put in the last block all monomials of degree bounded by t — d n+ i. 
The rows correspond to TZi for 1 < % < n. We also index them by graded 
descending order: first all monomials of degree t — di in TZi for 1 < i < n, 
then all monomials of degree t — di — 1 in TZi, 1 < i < n , an d so on up to all 
monomials of degree t — di — d n+ i + 1 in 7£j, 1 < i < n. In the last block we 
put all monomials of degree bounded by t — di — d n+ i in TZi, 1 < i < n. 

With this ordering M' has a block structure: 

M t * * * 
* * 

M = 

E 

where the square matrix Mj corresponds to the coefficients of the terms of 
degree j of x. a fi, where \a\ = j — d iy that is, the coefficients of x a /j except 
those corresponding to terms in Tj. 



(14) 
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Hence Mj is the Macaulay-Chardin matrix associated to the j-subresultant 
of f ± , . . . , f n with respect to Tj ([5]) and it turns out that 

|M,|=^A rj , 

where £j is the extraneous factor associated to this construction, that we recall 
only depends on j and not on the set Tj. 

But it turns out that the extraneous factor £(t) has a block structure similar 
to (14) (see [24,4,8]). We have, with our notation: 

£{t) = \E\ fl (15) 

j=t— dn+l+l 

(see [24, Th. 6]). This concludes the proof of the Theorem. □ 

Remark 3.3 The reason why we cannot allow Tj to be any subset of mono- 
mials of degree j for j < t — d n+ i + 1 is the factorization formula on the 
right hand side of (15), where the £/s involved in the product are only those 
corresponding to j satisfying t — d n+ i + 1 < j < t. This is not just a technical 
obstruction. If we could pick any Tj for every j, then setting t := p + d n+ i, the 
Poisson formula for the resultant Res(/f, . . . , f^ +1 ) would read as follows 

ef 1 • • • d 1 



SI ' ' ' Sd 



Res(f f n ) d ^ /n + i(0, 

T CeV^n (f U ...,fn) 

which is obviously false in general since the fraction does not cancel unless 
T = 1Z n+ i, i.e. Tj is defined as in (10). 

Like in the univariate case, we illustrate Theorem 3.2 with a specific example. 
Example 3.4 Let n = 2, d\ = d 2 = d 3 = 2 and t = t* = 2. 
Here k = #{xi,x 2 , Xix 2 } = 3, TZi = TZ 2 = 7^ 3 = {1} and T = {1, xi, x 2 , Xix 2 }. 
We fix the ordered monomial basis (1, xi, x 2 , x±x 2 , x\ , x|) of -^[x] 2 and 

fi — do + diXi + a 2 x 2 + 03X1X2 + + 05X3 

f2 = b + 61X1 + 6 2 X 2 + &3^1^2 + ^4^1 + &5^2 
h = C + C1X1 +C 2 X 2 + C3X1X2 + C 4 x\ + C5X2. 
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Then 



M h 




a ai a 2 a 3 a 4 a 5 


M h 




b 61 6 2 63 64 b 5 


M f3 




Co Ci c 2 c 3 c 4 c 5 



We choose 5 := {xi ) xix 2) x\}. Then 

A 5 = c (a 2 b 5 - a 5 b 2 ) - c 2 (a b 5 - a 5 6 ) + c 5 (a & 2 - a 2 b ). 

On the other hand, if V R (f 1 J 2 ) = {6,6,6,6} with 6 = (61,62) for 1 < 
j < 4, then 

61 61 61 61 
6162 6162 6162 6162 

tl f2 t2 tl 

Sll S21 S31 S41 

/ 3 (6) / 3 (6) / 3 (6) / 3 (6) 



o s = 



Therefore, if we set V for the generalized Vandermonde matrix on 6, • • • ,6 
corresponding to the sequence of monomials 1, x±, x 2 , Xix 2 , xf, x 2 ,, i.e. 



V 



1 


1 


1 


1 


11 


61 


61 


61 


'12 


62 


62 


62 



6162 6162 6162 6162 



6 2 1 

f 2 

S12 



£fl 

f 2 
S22 



f 2 
S 3 1 

f 2 
S 3 21 



6 2 i 
£42 



e K , 



and Vij, < i < j < 5, for the square submatrix obtained from V deleting 
the i-th and j-th rows (we adopt the convention of numbering the rows from 
to 5 like the coefficients of the /i's), we conclude that 

\O s \ = -c |V2 )5 | + c 2 |Vb, 5 | + c 5 I Vb, 2 1 • 

Also, with this notation is the Vandermonde matrix corresponding to T. 

Now, since the only non-trivial homogeneous subresultant A7- in (13) is for 
T 2 = {x\X 2 }, and is equal to 

Ar 2 = a 4 b 5 - a 5 6 4 , 



Theorem 3.2 states that 
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Co(a 2 b 5 - a 5 b 2 ) - c 2 (a b 5 - a 5 b ) + c 5 
= ±(a 4 b 5 - a 5 b 4 ) ( -c -|^ 



^0,5 ^ 
C 2777 ~ T + c 5 



■i - a 2 b ) 

V o\\ 



1^4,51 " | ^4,5 1 | ^4,5 1 / 

Indeed, we show below that this equality holds since for any i < j and k < I: 



(-1)*+* 



Wi 



■';J\ 



CLibj — CLjbi 



(-1) 



k+l 



\V k 



a k bi - aib k 



(16) 



If for < i, j < 5, we set Tjj G i^ 4xb a 4-identity matrix with added columns 
for column i and column j, and G K 6x2 the matrix with 4 null rows and 
the identity matrix plugged in rows % and j, we observe that 






* 





a k ai 
b k k 



since f\{ij) = /2(£j) = 0, 1 < j < 4. Thus, taking determinants on both sides, 

( _ 1)5 _ j+ 4- i(Q . fe . _ ^ . ( _ 1)fc+ ,_i| yfc>j| = . {akbi _ aihk) ^ 

and we obtain (16). 

Applying this to our case, we conclude that here 

|Os| 



n a 



\j=t-d n+ ±+l 



Vr 



□ 



\ r . = < 



Next, we recover Theorem 2.2 in the univariate case: 

Observation 3.5 For n — 1, by setting /i := g and / 2 := /, as f\ = bd 2 x d2 , 
it turns out that 

b d2 if j > d 2 
1 if j < d 2 . 

So, if t > d 2 , then X[) =t _ dl+l A Tj = b d .^ d2+x . lit < d 2 , the product of subre- 
sultants equals 1. □ 

In the particular case t — p + <i n +i, Theorem 3.2 gives a new proof for the 
Poisson product formula for the multivariate resultant (see [6]): 
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Corollary 3.6 (Poisson product formula) 

Res(fl . . . , O = ±Res(J 1 , . . . II WO- 

Proof. We apply Remark 3.1 (2) for t := p + d n+ i to Theorem 3.2. We 
observe that by the same remark, for j > p, i.e. for j > t — d n + 1, A^. = 

Res(/ 1 ,...,/J. We conclude that O s equals (ll^v 7r (h,...,f n ) /n+i(0) times 
the generalized Vandermonde matrix whose determinant equals Vr- n 



We end this paper by giving the multivariate version of Corollary 2.7, i.e. a 
discrete Wronskian type expression for the generalized subresultant polynomial: 



»W==EV- ( 17 ) 
j"=o 

defined for a fixed t e N and k := 7~td 1 ...d n+1 (t), under the usual notations, 
S := {x^ , < j < k} C fT[x]t and 5,- := S \ {x»}. 

It turns out that s(x) belongs to the ideal generated by the /j's (see [5]), and 
the following result can be proved mutatis mutandis the proof of Corollary 
2.6. 

Corollary 3.7 Let /i, . . . , / n +i G if [x] and s(x) fre £/ie generalized subresul- 
tant polynomial defined in (17). Then, we have 



S (x) = iVr- 1 | fl A Tj j 



X 



70 



X 



71 



£7o 



d 

71 



x 7fc 

£7n+l(£l) 
ei7n+l(6) 



t7fc 

£d/n+l(£d) 
£d7n+l(£d) 



□ 



Remark 3.8 If gcd(<S) G 5, then one can reduce the previous determinant, 
as in Corollary 2.7. 
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